b0=200
b1=100
set.seed(1)
x=seq(0,41)
fd=b0+b1*x
lines(fd)
comp=c(.07,.04, .11,.07,.04,.11,.11,.45)
N.init=array(0,dim=c(41,8,3))

for(j in 1:3){
for(i in 1:41){
	N.init[i,,j]=fd[i]*comp*runif(8,.5,1.5)
	}
}
N.init[1,,]=NA
#par(mfrow=c(1,1))
#matplot(N.init[,,3])
#save(N.init,file="/Users/Tom/Documents/YNP_Bison_Model/A_final_model_runs_for_revised_manuscript/Initial_state.Rdata")